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Abstract 

Increasing evidence of the effects of changing climate on physical ocean condi- 
tions and long-term changes in fish populations adds to the need to understand 
the effects of stochastic forcing on marine populations. Cohort resonance is 
of particular interest because it involves selective sensitivity to specific time 
scales of environmental variability, including that of mean age of reproduction, 
and, more importantly, very low frequencies (i.e., trends). We present an age- 
structured model for two Pacific salmon species with environmental variability in 
survival rate and in individual growth rate, hence spawning age distribution. We 
use computed frequency response curves and analysis of the linearized dynam- 
ics to obtain two main results. First, the frequency response of the population 
is affected by the life history stage at which variability affects the population; 
varying growth rate tends to excite periodic resonance in age structure, while 
varying survival tends to excite low-frequency fluctuation with more effect on 
total population size. Second, decreasing adult survival strengthens the cohort 
resonance effect at all frequencies, a finding that addresses the question of how 
fishing and climate change will interact. 
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i 1. Introduction 



In population theory, interest is increasing in the complex ways in which 
age-structured, density-dependent populations respond selectively to different 
time scales (or frequencies) of variability in the environment. These questions 
dovetail with the increasing practical need to understand how populations will 
respond to potential changes in life history parameters and the time scales of en- 
vironmental variability due to climate change and increasing pressure on natural 
resources such as fisheries. There is a growing awareness that model population 
responses can appear to filter certain frequencies in environmental variability 
(jGreenman and Bentonl 120051 ). Furthermore, these effects can serve to amplify 
variability by exciting modes of pop ulation behavior that without e nvironmental 



2 variability would be locally stable (jGreenman and Benton 



2003). One exam- 



ple of such behavior, cohort resonance, is closely identified with life history 
characteristics in that it features cycl es of period equal to the dominant age of 



5 reproduction. (jBiOrnstad et al 



1999) 



There is particular interest in the dynamic responses of the marine fish tar- 
geted by global fisheries as they are subject to a dramatic artificial change in a 
life history parameter, adult survival, and there are indications that the domi- 
nant time scales of environmental variables such as ENSO that affect fi sh popula- 



tions may change with a changing climate (jTimmermann et all 119991) . Empiri 



2i cal observations and mo dels indicate that variability in fish popula tions increases 



22 with increased fishing (jHsieh et al 



2006 



Anderson et al 



2008). The cohort 



resonance phenomenon in fish populations presents particular problems for cli 
mate change because it enhances sensitivity to very slow signa ls (trends) in ad- 



25 dition to those at the period of dominant age of reproduction (jBiOrnstad et al 



20041) . Observations of such effects could be confounded with potential slow 



changes in the environment, and thus make it difficult to differentiate between 
them. 
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29 Within the general concer n for the combin ed effects of fishing and climate 

30 change on marine ecosystems (jPerrv et all l2010f) , there is a particular interest 

31 in the effects of the marine environment and fishing on population dynamics of 

32 Pacific salmon on annual to decadal time scales. Analyses of the influences of the 

33 ocean environment on Pacific salmon, whether statistical exam ination of covari- 

34 abilit y between population and environmental variables (e.g., 



Logerwell et al 



20031) or estimation of survival s to specific sizes and ages through ana l ysis o f 



coded wire tagging data (e.g., 



Coronado and Hilborn 



1998: 



Teo et al 



2009), 



commonly assume the variable ocean environment influences survival during 
the early ocean phase of this anadromous genus. In fewer cases environmental 
variab ility in the age of m aturation has been the dependent population variable 



(e.g., iPyper et al 



1999). 



41 These empirical findings raise questions regarding the relative roles of ran- 

42 dom variability in survival at various ages, and in the age distribution of repro- 

43 duction, in salmon population dynamics. Random survival is typically assumed 

44 to influence abundance directly, but the effects of varying age of spawning are 

45 not as clear, and there is a need to understand differences in population response 

46 to these two sources of random variability. 

47 Here we examine and compare the effects of environmental variability in 

48 survival and age of spawning on the magnitude and time scales of population 

49 variability. We are engaged in a study of the influence of the oceanographic 

50 environment on Pacific salmon, the GLOBEC North East Pacific program, 

51 part of the US Climate Change program. Earlier retrospective analyses in- 

52 dicated differences in the responses of the two California Current congeners, 

53 chinook salmon, Oncorhynchus tshawytscha, and coho salmon, O. kisutch, and 

54 hence we are interested in th e population dynamic differences between them 



55 ( Botsford and Lawrence , 



20021) . For instance, coho salmon collapsed synchronously 



56 along the coast between 1980 and 2000, and coho appear to have more high fre- 

57 quency variability than chinook in their catch reco rd. These species diff er in 



their age distribution of spawning and in other ways (jBotsford et al 



2005). We 



have explored some of the differences in probabilities of extinction of populations 
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with these spawning age distributions in re sponse to time-varying marine sur 



vival (Hill et al 



2003 



Botsford et al 



2005). In salmon, differences in maturity 



sched ule are due to differences in size dist ribution (and therefore prior growth 
rate) (jYound . 11999 [_ Vollestad et al. L 20041) . Similar effects occur in other non- 



64 salmonid species ([Day and Rowe 



2002). To avoid confusion, we note that the 



mechanisms studied here are neither (1) the indirect effect on survival of varying 
develo pment rate due to co nsequent variation in time spent at higher mortality 

(e.g., 



Moloney et al 



1994I ). nor (2) the interaction between age structure and 



over-compensatory density-depend ent recruitment that causes c ycles of period 
roughly twice the generation time ([Ricker , 



1954 



Botsford, 



1997). 



In the models used here oscillations with period 2T do not arise (where T is 
generation time or age of dominant effect on recruitment), primarily because of 
the form of stock-recruitment function we use, but the cyclic mode with period 



Myers et al 



T tur ns out to play a central role in our investigation. Recently 
(|1998f ) explored stochastic forcing of a cyclic mode of variability of period T as 
a potential cause of observed cycles in sockeye salmon ( Oncorhynchus nerka) in 
British Columbia, Canada. This mode is similar to the "echo effect" associated 



with linear age-structured models of semelparous species (|Svkesl . 119691 ). Myers et 



al. noted that while this mode would not be the dominant mode, it could appear 
clearly in solutions obtained through forcing by randomly varying survival. This 
was essentially cohort resonance, though they did not use that term. 

The time scales (or frequency content) of the population response also de- 
pend on the mode of observation , e.g., whether t he data in a time ser ies are 



recruitment, abundance or catch ( Botsfor' 



rd, 



1986 



Anderson et al 



2008). The 



84 nature of a catch time series depends on fishing mo rtality rate, with higher fish 



85 ing ra tes leading to higher frequencies of variability ([Botsford 



1986 



Hsiehetal 



2006T ). Differe nces between the s pectr a of recruitment and abundance have been 



87 illustrated by 



Biornstad et al 



([20041 ) for Atlantic cod {Gadus morhua). The 
most commonly available measurement of salmon populations is annual catch, 
either in numbers or biomass, and occasionally the age distribution of catch 
is determined. Frequently catch can be assumed to be individuals who if not 
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91 caught would be spawning that year because they are near the mouth of (or in) 

92 the spawning river. In some streams, spawning escapement is also estimated, 

93 and sometimes the age composition of spawners is also estimated. Catch and 

94 escapement can be summed to obtain the total abundance of potential spawners 

95 in a year. Total population abundance at time t cannot be observed directly, 

96 but estimates can be obtained through cohort reconstruction from several years 

97 of age-specific data from catch and escapement. 

98 Here we explore several aspects of the spectral response of marine fish popu- 

99 lation dynamics, using parameter values representative of two species of Pacific 

100 salmon over a range of survivals as examples. Our practical interests are in 

101 the combined effects of fishing and climate change, so we focus on effects of 

102 long-term changes in survival on the population response to various time scales 

103 of environmental variability. Note, however, that long-term declines in survival 

104 can also be caused by the climate, independent of the fishery, as occurred i n 



2005). 



105 California Current coho salmon between 1980 and 2000 ([Botsford et al 

106 We also explore the differences in responses to environmental forcing at differ- 

107 ent points in the life history (i.e. development rate vs. survival), and spectral 
los differences between various kinds of observations (e.g., observations of recruit- 
109 ment vs. observations of total abundance). We present results in terms of the 
no changing spectral responses, and relate the spectral response to the eigenvectors 
in and eigenspace structure of the linearized model. 

n2 2. Model Formulation 

in To investigate these questions, we introduce a density-dependent, stochastic 

n4 age-structured model which we can use to compare the effects of temporal vari- 

n5 ability in the timing of reproduction and annual survivorship of salmon. Let n 

no be the maximum spawning age in years, and let x(t) = (x\(t), . . . ,x n (t)) T be 
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the age-structured population vector at year t. Then our model is 



x(t) = F(x(t -!),*) = 



R(P(t)) 
s(t)xi(t - 1) 
a(t) x 2 (t-l) 



\ 



(1) 



s(t) x n ~ 3 {t - 1) 
(1 - 5 e {t)) s{t) x n - 2 (t-l) 
\ 5i(t)s(t)x n -i(t-l) J 
where P(t) = (S e (t) x n -2(t — l) + (l — 6i(t)) x n -i(t—l) + x n (t— 1)) is the number 
of fish returning to spawn in year t, and 

<xP(i) 



R(P(t)) 



1 + pP(t) 



(2) 



is a Beverton-Holt density dependent recruitment term (IBeverton and Hold . 
19571 ). Thus the recruitment R(P(t)) represents the number of outmigrants 
(smolts) leaving a spawning stream in year t, resulting from egg production by 
individuals of various spawning ages. We have assumed that migration from 
freshwater to the marine environment occurs in the first year, and that there 
is no difference in fecundity between the different ages of spawning. Parame- 
ters a and j3 characterize the density-dependent reproductive phase: a is the 
density-independent per-capita growth rate when the population is very small, 
and ^ is the maximum total number of offspring in the population. Timing 
of reproduction is controlled by the remaining two parameters, 8 e (t) and Si(t): 
most individuals spawn at age n — 1; a proportion 6 e (t) of those surviving to age 
n— 2 spawn at that age in year t — 1; and a proportion 8i(t) of age n — 1 fish in 
year t postpone spawning until year t+1, when they are age n. Annual survival 
s(t), as written here, affects all age classes in year t. However, we also explore 
the possibility that the dominant variability in ocean survival occurs during the 
period immediately following ocean entry. This is a challenging stage for juve- 
nile salmon, as they are completing the transition from freshwater to the marine 
environment, and they are dependent on the highly variable food production in 
the coastal ocean. 
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139 We are interested in understanding the time scales of variation in the solu- 

140 tions to the model (equation[I} in response to fluctuations in annual survivorship 

141 and timing of spawning. We do this by studying three cases: the case of fluctuat- 

142 ing survival s(t) at all ages, the case in which survival varies but only at the age 

143 of entry into the ocean, and the case in which fluctuating mean age of spawning 

144 a(t) produces yearly changes in S e (t) and Si(t). In all cases, the environmental 

145 fluctuation is a Gaussian white-noise signal, € R with E(£(t)) = 0. To 
«6 model fluctuating survival we use s(t) = s° + where s° is the unperturbed, 

147 constant survival. In the simulation the distribution of s(t) is truncated to en- 

148 sure that < s(t) < 1. To model fluctuating age of spawning, we suppose that 

149 individuals' ages of spawning are chosen from a normal distribution whose mean 

wo is the central age of spawning a(t) in year t (Fig. [TJ, specifically 

I f a-ait) y 
p(a-a{t)) = ^=^e I ^ J . (3) 

151 We approximate early and late spawning by having all early spawners spawn at 

152 age n — 2 and all late spawners spawn at age n, so that 

/n— 1.5 
p(a-a(t))da (4) 
-OO 

153 and 

Si(t) = / p(a-a(t-l))da. (5) 

Jn-0.5 

154 From the point of view of a cohort, the portions of that cohort that will spawn 

155 early and late, i.e., S e and Si, are set by the value of a(t) in the year that that 

156 cohort transitions from age n — 2 to age n — 1. That is, 5 e (t) is a function of the 

157 central spawning age in year t, whereas 6i(t), which has its effect a year after it 

158 is determined, is a function of the central spawning age in year t—1. We model 

159 the fluctuating central age of spawning as a(t) = a° +£(£). 

160 To compare the dynamics of the two salmon species along the west coast of 

161 the contiguous U.S., we have used parameter values that approximate known 

162 or likely values for the populations of coho and chinook salmon. Survival 

163 rates in the ocean are reasonably well known, and we consider three cases: 

164 a "typical" value of 0.85 yr 



Bradfor 



1995h . a "small" value of 0.28 yr~\ 
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n-3 n-2 n-1 n n+1 

Age 



Figure 1: An example illustrating how the fraction of early and late spawners is modeled. The 
central spawning age a(t) varies from year to year. Here we picture a normal distribution of 
potential spawning age with mean a(t) = n — 1.1 and a = 0.8. The parameters representing 
the proportions of early spawners (5 e ) in year t and late spawners (5j ) in year t + 1 are the 
integrals of this function over (-co, n — 1.5] and [n — 0.5, oo), respectively. 



165 and a "very small" value of 0.2 yr _1 , estimated from more recent observa- 

166 tions of coho salmon. Since the "very small" value is not sustainable in the 

167 chinook model, that case is not considered. The distributions of values of a 
tea and ft ( equation [2t have been es timated for coho salmon, but not for chinook 



salmon ( Barrowman et al 



2003), and we use the modes from those distribu- 



170 tions (a — 60 and f3 — 0.00017). In chinook salmon populations in the California 

171 Current, individuals typically spawn primarily at a single dominant age, with 

172 less spaw ning at ad j acent ages. The dominant age of spawning increases with 



20031 ). Here we chose dominant spawning at age 4, i.e. 



173 latitude (jHill et al. 

174 n = 5 and a° = 4. Precocious spawning in chinook salmon range s from 0- 



1991). We use 



175 10 percent in males and from 0-3 percent in females ([Healev . 

176 a — 0.4, which makes 5 e = 8i w 0.1. Coho salmon in the California Current 

177 spawn predominantly at age 3 with substantial precocious spawning at age 2 

178 and minimal spawning at age 4. Precocious spawning in coho salmon consists 

179 almost completely of males, is more variable than i n chinook salmon, and can be 
iso as high as 30 percent in natural wild populations (|Sandercocld . 119911 ). It is not 
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181 possible to observe the effects of precocious spawning by males on reproduction 

182 directly, but for coho salmon, a modification of genetic methods for estimating 

183 effective population size indicated the effective pro portion of 2-year-olds to be 



2002). For 



184 35 percent in two naturally spawning populations (jDoornik et al. 

185 coho salmon we chose n = 4, a = 2.75 and a = 0.2, which makes S e w 0.1 and 

186 Si « 0.0001. The standard deviation of the stochastic variation, erg, is 0.085 for 
variation in survival, and 0.2 for variation in mean spawning age. 



Table 1: Model parameters and unperturbed values. 



Parameter 


Description 


Coho Value 


Chinook Value 


s° 


annual survival rate, all ages 


0.85, 0.28, 0.2 


0.85, 0.28 


Se 


proportion of age n— 2 fish spawn- 
ing 


0.1056 


0.1056 


Si 


proportion of age n — 1 fish delay- 
ing spawning to age n 


8.84 x 10" 5 


0.1056 


a 


slope at origin of Beverton-Holt 
stock-recruitment curve 


60 


60 


P 


saturation parameter of 
Beverton-Holt stock-recruitment 
curve 


1.7 x 10~ 4 


1.7 x 10~ 4 


a° 


central age of spawning 


2.75 


4 


(T 


standard deviation of spawning 
age distribution 


0.2 


0.4 


CTE 


standard deviation of environ- 
mental forcing signal 


0.085 (survival) 
0.2 (spawning age) 


0.085 (survival) 
0.2 (spawning age) 



las 3. Model Analysis 

189 In each of the model cases, the age vector x converges to a neighborhood of a 

190 single point, and fluctuates in that neighborhood in response to the fluctuating 

191 life-history parameters. To find the center point, we consider the deterministic 

192 system defined by equation [1] with s, 5 e , and Si all fixed at their unperturbed 

193 values. This system has a unique positive fixed point that is globally attracting 
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194 for positive trajectories: 



s 

<,2 



„n-3 



(6) 



where 1/c 



n— 3 



(1 - 5 e )s n - 2 



+ (1 - <5 e )(l - (50s"- 2 + (1 - 8 e )Sis n - 1 is the amount of 



spawning in the lifetime of the average smolt. 

As we will see, we can understand a great deal about the population dynam- 
ics by looking at the linearized dynamics at this fixed point of the deterministic 
dynamics. We linearize the effect of the stochastic noise term as well as the 
age-structured population vector x(t): 

L 

jftt) « Jjfti-l) + -0. (7) 

1=0 

where y(t) = x{t) — x{t) is the deviation from the fixed point, J is the Jacobian 
matrix of the dynamics at the fixed point, and linearizing in the noise term gives 
a sequence of "forcing vectors" H that summarizes how the noise enters into 
the dynamics with various time lags (see the Appendix for details). 
At the fixed point, the Jacobian matrix is 

/ ••• 5 e c - - \ 

e a v L > a. a 

s 



J = 



(1 - 6 e )s 



(8) 



5is J 



206 where all entries left blank are zero. 
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A linear system like this one decomposes naturally into independent subsys- 
tems located in linearly independent one- and two-dimensional s ubspaces, each 



209 chara cterized by a real eigenvalue or a complex pair of eigenvalues (jHirsch and Smale 



1974T ). The right eigenvectors of the Jacobian matrix are the basis vectors for 



211 each of these subspaces of the dynamics. The complex conjugate eigenvalues 

212 of the Jacobian predict the resonant frequencies of the population's response to 

213 noise. In the time-varying-survival cases, we only need a single forcing vector 
2« H, because the state of the environment in year t, only affects survival in 

215 year t; but in the varying-age-structure case, because conditions in year t affect 

216 the number of early returns S e in year t and the number of late returns 6i in 

217 year t + 1, we have to include two forcing vectors to describe the effects with 

218 and without one year's time lag. 

219 In the case of fluctuating survival at all ages, the forcing is captured by the 

220 vector H s (0), with 



/ 



o 

X3 



%n—3 
(1 - S e )x n - 2 



\ ( 







(s°)"- 4 

(i-s e )( s °r- 3 
\ (i-5 e )^( S °)"- 2 j 



(9) 



221 We note that when survival is forced additively, as we have done here, the 

222 deterministic system described above is not exactly the mean of the stochastic 

223 system. A stochastic system with s(t) = s°e^^ would have the deterministic 

224 system as its mean, but in the limit of small noise, these representations have 

225 identically-shaped frequency responses that merely differ by a factor of s. We 

226 chose additive noise for convenience. 

227 For forcing of survival at ocean entry, only survival to age 2 is subject to 
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228 fluctuation, so that 



/ n \ 



a — c 




(10) 



V / 

229 For time- varying ages of maturation, we have lag effects from early spawn- 

230 ing 

= = (11) 



(12) 



^ (0) - \m) - Kdfemj 

231 and lag 1 effects from late spawning 

( dF t \ fdF dSi 



Ha(l) 

Since = 1 we have 



\dZ(t-i)J \d6id£(t-i) 



dS e (t) 



-I 



n-1.5 



-p'(a - a(t j) da = -p(n - 1.5 - a{t)) (13) 



and 



dSi jt) 



— = / -p'(a - a(t- l))do = f>(n - 0.5 - a(t - 1)). 

1) Jn-0.5 



234 The vector derivatives we need are 



a — c 



and 



m 

dS l 





-s"- 2 
V J 

-{l-8°)s n - 2 ^ \ 




a — c 



\ (1 - 5°).s«- 1 / 



(14) 



(15) 



(16) 
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236 where 6° refers to the unperturbed value of 5 e . Substituting (TT3"|) and (fT5|) into 

237 (fTTj) yields the forcing at lag 0, 

( -p(n- 1.5- a°)s n ~ 3 ^ \ 




i?a(0) = 



V 





p(n - 1.5 - a°)s n ~ 2 




a — c 



(17) 



/ 



238 Substituting (|14p and (jTHJ) into (fT2)) gives the forcing at lag 1, 

( -p(n-0.5-a°)(l-5°)s n - 2 ^ \ 




^a(l) 



(18) 



\ p(n-0.5-a°)(l-(5°)s"- 1 / 



239 4. Results 

240 ^.L Dynamics near the fixed point 

241 Each of these models without random fluctuation has a unique positive fixed 

242 point (equation ©)• The condition for a fixed point can be illustrated in terms 

243 of the stock-recruitment function and a straight line through the origin with 

244 slope equal to the inverse of mea n lifetime egg production, as shown in Fig. [5] 



245 (jSissenwine and Shepherd 



19871) . The fixed point lies at the intersection of 



246 these, allowing clear interpretation of the effects of long-term changes in sur- 

247 vival. As mean ocean survival declines, equilibrium egg production declines, 

248 and eventually recruitment declines. For the very low survival rate for chinook 

249 salmon, the fixed point is zero indicating survival is not adequate for persis- 

250 tence. Because the parameters for the stock-recruitment relationship and the 

251 estimated survivals are for coho salmon, this has no implications for real chinook 

252 salmon populations. 
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Figure 2: Graphical interpretation of the equilibrium recruitment as the intersection of the 
smolt-adult curve and a line through the origin with slope l/(lifetime reproduction) for (a) 
chinook and (b) coho salmon. 
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-1 1 

Re(J-) 

Figure 3: Eigenvalues of linearization matrix J in chinook and coho models for typical survival: 
(a) chinook, s = 0.85 yr" 1 , (b) coho, s = 0.85 yr — x ; small survival: (c) chinook, s = 0.28 yr — x , 
(d) coho, s = 0.28 yr — 1 ; and very small survival: (e) coho, s = 0.2 yr — 1 . 

253 Plots of the eigenvalues of the linearization at the fixed point of each of 

254 these models indicate the modes of variability, and how they change with sur- 

255 vival (Figure [3]). Note that the chinook model has 5 eigenvalues, while the coho 

256 model has 4 eigenvalues because those are the maximum ages in each model 

257 population. In both the chinook and coho models, the linearized population 

258 dynamics respond most strongly to forcing in the one-dimensional mode corre- 

259 sponding to the positive real eigenvalue, Ai (the dominant mode of variability). 

260 Motion in this mode occurs without oscillation, that is, at low frequencies. The 

261 next strongest mode is a resonance at period n — 1 or a little less (determined by 

262 the eigenvalues A2 and A3, whose angle from the positive axis is approximately 

263 ±-^fr), indicating oscillations with a period equal to one generation time. In the 

264 chinook case (Figures 3a, c) there is also a strong resonance at period two, corre- 

265 sponding to a negative eigenvalue, and in all cases there is a weakly resonating 

266 negative eigenvalue as well. 

267 In both the chinook and the coho cases, the effect of decreasing the mean 

268 survival is to move all eigenvalues outward toward the unit circle, increasing the 
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269 return time of all subsystems of interest and thereby increasing the expected 

270 magnitude of the cumulative response to variation over time. This effect is 

271 especially strong in the chinook model for "small" s (Fig. [3};). This occurs 

272 because the rate of change of recruitment with adult stock increases substantially 

273 as the survival declines (Fig. [2j. Therefore, low equilibrium abundances in 

274 Figure [5] correspond to large magnitude eigenvalues (and, thus, large responses 

275 to variation)in Figure[31 

276 The spread in ages of spawning has a relatively small specific influence on 

277 the locations of the eigenvalues. When a — 0, all fish spawn at age n — 1, 

278 and the eigenvalues are roots of a simple characteristic polynomial, P(X) = 

279 A(A ,l_1 — s™~ 2 ^-;). There is one zero root and the others are equal in size 

280 and evenly spaced around zero (not shown). As a increases, the dominant, 

281 positive real eigenvalue moves outward on the negative real axis, the others 

282 move inward slightly and the complex eigenvalues rotate slightly, their polar 

283 angles becoming slightly smaller in the chinook case and slightly larger in the 

284 coho case. In all cases other than a = 0, the positive eigenvalue is largest in 

285 magnitude (consistent with the Perron- Frobenius theorem, since the entries of J 

286 are nonnegative) . For our purposes, we conclude that the spread in spawning age 

287 has very little impact on the model dynamics (too slight to justify an illustration, 

288 in fact), and so our models do not indicate that it is an important difference 

289 between chinook and coho population dynamics. 

290 1^.2. Mechanism of environmental forcing 

291 To illustrate the difference between population responses to different mech- 

292 anisms of environmental forcing we compare the cases with (1) variability in 

293 survival in each ocean year, (2) variability in survival at the age of ocean entry 

294 only and (3) variability in the spawning age distribution. Random variabil- 

295 ity in survival appears to preferentially excite the geometrically decaying mode, 

296 while variability in the spawning age preferentially excites the cyclic mode. This 

297 conclusion holds for all five of our model cases: the resonant frequency is about 

298 0.33 yr _1 for the coho salmon, whose dominant age of spawning is three years, 
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299 and about 0.25 yr -1 for the chinook salmon, whose dominant age of spawning is 

300 four years (not shown here, see next example). The one exception is that for the 

301 chinook model with s° = 0.28 yr -1 , there is a strong low-frequency component 

302 with variation in spawning age as well as the period-4 component. 

303 We illustrate this general result with the case for coho salmon with s° — 

304 0.28 yr -1 by comparing the frequency responses to survival varying at all ages, 

305 for survival varying only at age 1, the presumed age of ocean entry, and for vari- 

306 ability in spawning age, then presenting examples of time series for each case. 

307 Both frequency responses to environmental variability in survival (Fig. |4^,) de- 

308 cline from very low frequency, leveling off slightly at a frequency just below that 

309 corresponding to period 3 (the dominant age of spawning for coho salmon) , then 

310 decline for higher frequency. The simulations indicate up to approximately 20 
3n percent greater variability than the analytical model in both survival cases. The 

312 case with variable ocean survival at all ages (Fig. \%jp) is skewed more toward 

313 variability at low frequencies. The frequency response to environmental variabil- 

314 ity in the spawning age distribution increases from low frequency to a resonance 

315 at a frequency slightly greater than that corresponding to period three, the 

316 dominant age of spawning, consistent with cohort resonance. 

317 The time series for each case (Fig. [5j indicate a visually discernable differ- 

318 ence between populations driven by environmental variability in survival and 

319 environmental variability in spawning age distribution. The population with 

320 environmental forcing of spawning age (Fig. [St) has clear indications of cohort 

321 resonant behavior near period 2 yr to 4 yr, while the populations with environ- 

322 mental forcing of survival tend toward 5-10 yr fluctuations with little variability. 

323 The difference in the magnitude of variability between the series with varying 

324 survival is as expected from the difference in area under the two thin lines in 

325 Fig. Hit. 

326 ^.3. Population observation 

327 To determine the effects of the type of observation on the time scales of 

328 variability we compared the results of observing recruitment to the results of 
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Figure 4: Magnitude of frequency response of recruitment for the linearized model (thin 
lines) and nonlinear simulations (bordered lines) for environmental forcing of (a) early ocean 
survival rate (solid lines) and ocean survival at all ages (dashed lines), and (b) spawning 
age distribution for coho salmon with s° = 0.28 yr — x . See the Appendix (equation lA.lOl l 
for derivation of the linearized frequency reponse. Nonlinear frequency response curves are 
estimated by calculating the average magnitude (the square root of mean power at each 
frequency) of the fast Fourier transforms of 1000 simulation time series of length 128, and 
scaling the result by \/T28 erg, to obtain the same units as the analytically-derived transfer 
function. 

329 observing total abundance for each case. Generally, they all exhibited rcla- 

330 tively more 2-to-5-year oscillation in the recruitment and more low-frequency 

331 fluctuation in the total abundance. 

332 To illustrate the difference between different types of population observa- 

333 tions, we show the results of observing recruitment with the results of observing 

334 total abundance, for chinook salmon with typical survival (s = 0.85 yr -1 ) and 

335 environmental variability in the spawning age distribution (Fig. [6]) . The spec- 

336 tral response of recruitment increases from low frequency to a peak at about 

337 0.25 yr -1 , the frequency expected for chinook salmon with dominant age of 

338 reproduction at 4 yr, then declines. The spectral response of total abundance, 

339 the sum of several cohorts, declines monotonically from low frequency, showing 

340 only a hint of the resonance present in recruitment. 

341 The time series of these two cases (Fig. [7]) reflect these characteristics. The 

342 time series from the population with observations of recruitment of a chinook 

343 salmon population with environmental forcing of spawning age appears to have 

344 a preponderance of variability on time scales of 2 yr to 5 yr, while the time 
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Figure 5: Time series of recruitment in simulations with environmental forcing of (a) early 
ocean survival (<te = 0.085), (b) ocean survival at all ages (erg = 0.085), and (b) spawning 
age distribution (ag = 0.2) for coho with s° = 0.28 yr -1 . Identical sequences of standard 
normal random variables were used as the forcing signal in order to allow a comparison of 
filtering by different demographic mechanisms. 




Figure 6: Magnitude of frequency response of (a) recruitment and (b) total abundance for 
the linearized model (black line) and nonlinear simulations (bordered line; 1000 simulations 
of length 128) for environmental forcing of spawning age distribution for chinook salmon with 
s = 0.85 yr" 1 . 
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Figure 7: Time series of (a) recruitment and (b) total abundance from a single simulation 
with environmental forcing of spawning age distribution (erg = 0.2) for chinook salmon with 
s = 0.85 yr" 1 . 

345 series from the observation of abundance appears to be a smoothed version of 

346 that. 

347 4-4- Long-term mean survival 

348 To illustrate the differences between populations operating at different levels 

349 of long-term survival we compared the coho salmon recruitment from a model 

350 with environmental variability in spawning age distribution with each of the 

351 three survival levels, typical, small and very small (Fig. [5]). The shift to lower 

352 constant survivals could result from fishing or a shift in climate. The frequency 

353 responses with these three survivals have similar shapes but differ substantially 

354 in scale. The resonance at period just less than 3 years is dominant because 

355 the variability is in spawning age distribution. Note the greater increase in 

356 sensitivity to low frequencies as the survival decreases to its lowest level. 
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Figure 8: Magnitude of frequency response of recruitment for the linearized model (thin lines) 
and nonlinear simulations (bordered lines; 1000 simulations of length 128 with = 0.2) for 
environmental forcing of spawning age distribution for coho with s = 0.85 yr — 1 (solid lines), 
0.28 yr" 1 (dashed lines), and 0.2 yr~ x (dash-dot lines). 



357 The time series of these three cases (Fig. [§]) appear to have similar frequency 

358 content, but different levels of variability as expected from Fig. [8] Importantly, 

359 they also underscore the fact that as survival declines variability increases as in 

360 Fig. [51 as the equilibrium recruitment declines. 

361 5. Discussion 

362 These analyses provide an understanding of how various population charac- 

363 teristics shape the response of salmon populations to environmental variability 

364 on various time scales. Salmon population variability does not simply follow 

365 variability in the environment as is commonly assumed, but rather the observed 

366 response is shaped by three factors: (a) the life history point of impact of 

367 the environment, (b) how the population is censused, and (c) pre-conditioning 

368 by long-term changes in survival. These considerations have important conse- 

369 quences for the management of salmon populations and the anticipation of the 

370 effects of large-scale environmental change. The dynamics of salmon responses 

371 to the environment are, of course, not unique, but rather are closely related to 



372 those of other higher t rophic level, age-structured species ([McCann et al 



Biornstad et al 



2004) 



2003 
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Figure 9: Time series of recruitment in simulations with environmental forcing of spawning 
age distribution (erg = 0.2) for coho with s = 0.85 yr — 1 (solid black), 0.28 yr — 1 (white with 
black border), and 0.2 yr -1 (dashed white on black). Identical sequences of standard normal 
random variables were used as the forcing signal in order to allow a comparison of filtering at 
different survival rates. 



374 The importance of cohort resonance to climate change is enhanced by the 

375 finding that its effects intensify with decreasing survival. There is concern over 

376 the effects of climate change on fisheries, and ho w management should change to 

377 mitigate those changes {e.g. 



Perry et al 



2010|) . The basic resonance mode with 



period T appears as the second and third eigenvalues in the linearized model, a 
conj ugate pair. It is well known in linear age-structured models, as the echo ef- 



fect (Svkc; 



1969 



Caswell . 



20011 ). and was identi fied earlier as a pot ential cause 



38i of observed cycles in some salmon populations ([Myers et al 



199S). However, 



382 increasing cohort resonance with fishing has not been mentioned as a factor in 

383 the variability in fished populations. Knowing that high levels of fishing and 

384 specific time scales of environmental variability can increase overall population 

385 variability will be valuable in formulating management policies to maintain sus- 

386 tainability. The additional variability increases the risk of population collapse 

387 in addition to that risk incurred by reducing mean abundance. In addition, 

388 knowing that fish populations could be more sensitive to some time scales than 

389 others will aid in the explanation of changes in fisheries' variability if the time 

390 scales of forcing change with climate. Potential intensification of cohort reso- 

391 nance through increased fishing also underscores the importance of Bj0rnstad's 
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392 (|2004r ) warning regarding the low-frequency effect of cohort resonance, sensi- 

393 tivity to very low-frequency environmental variability. Fishing increases the 

394 chances that populations, even when driven by white noise in the environment, 

395 could generate very slowly changing signals that could be mistaken for the ef- 

396 feet of a slowly changing climate. This hard-to-detect long-timescale change 

397 may present a third increase in threat when coupled with reduced abundance 

398 and increased variability due to fishing. Note that this increased sensitivity at 

399 low frequencies was seen in coho salmon (Fig. 8) and also in chinook salmon 

400 with s = 0.28 yr _1 (as noted at the end of section 4.2). 

401 Differences in the dynamic responses of populations to temporal variability 

402 at different points in the species's life history are not commonly considered. For 

403 salmon (and most other marine fish), while the effects of variability in survival 

404 are widely appreciated, the authors are unaware of any description of the ef- 

405 fects of variability in spawning age distribution on population dynamics, such 

406 as on cohort resonance. Here we found that the different mechanisms involving 

407 variability in survival and variability in spawning age distribution preferentially 

408 excited different fundamental modes of variability; hence these differences are 

409 clearly important. In addition to variability in spawning age exciting a com- 

410 pletely different mode than variability in survival, responses to varying survival 
4n at ocean entry and varying total survival also differed, with less sensitivity to 

412 low frequencies resulting from variability in early ocean survival than from vari- 

413 ability in total survival. This further motivates efforts to discover the ages and 
4« locations at which variability in salmon marine survival occurs. This may be 

415 related to the difference in po pulation persistence between these two survival 

416 mechanisms in earlier analysis (jBotsford et all 120051 ) , where population variabil- 

417 ity was greater when it occurred at the age of return to the river for spawning, 

418 rather than the age of ocean entry of smolts. In that case this difference occurs 

419 because of the Law of Large Numbers, and the fact that the logarithm of the 

420 number of spawners is the sum of several random survivals when variability is 

421 at the age of entry, but only one random survival when it is at the age of return. 

422 In our models, varying age of spawning causes much more oscillation at the 
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423 period of the generation length and less low-frequency fluctuation than varying 

424 survival does. This might be because changes in central spawning age have a 

425 double effect on recruitment, removing fish from one year's stock of potential 

426 spawners and adding them to the next or previous year's potential spawning 

427 stock. This perturbation has a direct effect on two cohorts' recruitment, which 

428 is echoed in subsequent generations when those cohorts spawn. Fluctuations 

429 in survival, on the other hand, whether they affect one cohort or many, cannot 

430 have this double effect, with the likely consequence that the period- T echo effect 

431 is not as extreme. 

432 The difference in frequency content between a c atch (or abunda nce) series 



Botsford 



19861 ). but there 



433 and a recruitment series was known previously (e., 

434 has been renewed interest in it. That difference basically follows from the Law 

435 of Large Numbers and the fact that when recruitment is the source of temporal 

436 variability, abundance will have a lower coefficient of variation because it is the 

437 sum of several recruitments, and hence can be expected to be a smoothed version 

438 of the recruitment signal, i.e., to emphasize lower frequencies. As the value of 

439 constant harvest rate increases, the number of cohorts summed diminishes, and 

440 hence higher frequencies are observed. 

441 The major difference between the population dynamics of coho salmon and 

442 chinook salmon revealed here is that they would be most sensitive to different 

443 time scales of environmental variability, as determined by the difference between 

444 the dominant age of spawning in each. Also, one explanation proposed for the 

445 fact that coho salmon declined near 1980, while chinook salmon did not, was 

446 that narrower width of the coho spawning age distribution. The effect of width 

447 of the spawning age distribution we determined here was not strong enough 

448 to support this hypothesis, leaving open the possibility that coho salmon were 

449 merely following a decline in ocean survival which did not decline for chinook 

450 salmon. These theoretical results are not definitive regarding specific stocks; 

451 rather they provide a context for further detailed investigations of the hundreds 

452 of coho and chinook salmon populations along the west coast of the U.S. 



As with other structured-population models (jNisbet and Gurneyi . ll982l : lBiOrnstad et al 
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2004; 



Greenman and Benton 



20031) . this model's response to environmental 

455 fluctuation is predicted well by the linearization at the deterministic model's 

456 fixed point. Like Greenman and Benton, we find amplification of noise as pa- 

457 rameters vary — in our case, environmental noise is amplified more by the 

458 population dynamics when survival is reduced. In our case, however, this am- 

459 plification is not associated with a nearby bifurcation, only with eigenvalues 

460 moving closer to the unit circle, a more general phenomenon. 

461 As discussed above, we also observe major differences in frequency responses 

462 among different kinds of environmental influences and between different popu- 

463 lation measurements. The analysis in our Appendix explains these differences 

464 in terms of different geometric relationships between the forcing vectors, the 

465 eigenspaces corresponding to the eigenvalues of the linearization, and the mea- 

466 surements. The projection of the forcing vectors H into each dynamical sub- 

467 space determines how strongly the environmental noise stimulates motion in 

468 that subsystem, and therefore how prominent resonance at that frequency is in 

469 the population dynamics. Similarly, different observations — total population 

470 and recruitment — include the different subsystems in different proportions, 

471 and so each mode of population behavior is more visible in some observations 

472 than others (see the Appendix for the mathematical treatment of these ideas) . 

473 This analysis, together with the diversity in frequency responses we see in our 

474 models, points to the importance of understanding the relationship between the 

475 eigenspace structure of the linearization and the relative importances of the 

476 resonant frequencies associated with each eigenspace. As an analytical tech- 

477 nique, this approach may be useful in multispecies ecological models as well as 

478 in structured population models. 

479 The results obtained here provide the means both to begin to explain current 

480 differences in responses to the environment by the same species at different 

481 locations, and to project differences in future responses on the basis of projected 

482 changes in time scales of variability in the environment. Current differences in 

483 responses by populations of the same species are typically presumed to imply a 

484 difference in environmental forcing, but they may be due to differences in life 
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485 histories or differences in preconditioning because they are fished at different 

486 intensities. Future changes in time scales of variability are expected on the 

487 basis of paleological records (e.g., past changes in the ti me scales of variabi lity 

488 of El Nino events as observed in corals and other media, 



Jones et al 



1999). 



2009) or 



489 predictions from global climate models (jTimmermann et al 
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609 Appendix A. General Mathematical Results 

610 The linearization of a stochastic map 



x(t) = F(x(t-l),^t),^(t-l),...,^t-L)), f(t)el", £(t)€R, (A.l) 



en where (without loss of generality) E(£(i)) = 0, is 



i7(t)« jtf(t-i) + 2#(0£(t-0 



(A.2) 



where x(t) = x + y(t), x is a fixed point of F(x, 0), J = (ff^-j |(=j ) ^ s * ne 
Jacobian matrix of F at x = x and £ = 0, H{1) = q^—^(x,0) is a vector 
expressing the dependence of F on the noise terms, and L is the maximum 
time lag at which stochastic perturbations affect F directly. Let Vi and Uj be 
the left and right eigenvectors of J, respectively, and A, its eigenvalues, so that 
ViJ = XiVi and Jili = In this paper we only consider matrices that have 

all distinct eigenvalues. 

We change to the natural coordinate system of J: Let 



U = (ui, . . .,u„), V 



I* \ 



( Ai 








\ 





A„ J 



(A.3) 



with J = UAV = Y,i Aifijtfi and UV = VU = I. Then for w = Vy, 

L 

1(t) = Vy(t) = AVy(t - 1) + £ 6{l)^t - I) 

1=0 
L 

= Aw(t - 1) + G(lK{t - I) 



w( 



(A.4) 



with G(l) = VH(l). Using the terms of this transformed vector w, the vector 
of deviations from equilibrium can be written as a sum of eigenvectors, y = 
J2i w$i- Since A is diagonal, the dynamics of each Wi{t) is uncoupled from the 
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624 others: 

L 

Wi (t) = \ iWi {t-l) + J29i(l)Z(t-l) (A.5) 
1=0 

625 where = ViH(l) is the iih entry of G(l). Thus, w(t) is the state of the 

626 linear system decomposed into its independent subsystems, and the vectors 
62? G(l) represent the stochastic forcing resolved into the decomposed coordinate 

628 system. 

629 Forcing and measurement 

630 The entries of the transformed vector G reveal how strongly the environmen- 

631 tal forcing acts on each subsystem of the linearized system. Different kinds of 

632 forcing (i.e. survival, age of maturation) are distributed differently among the 

633 different subsystems, characterized by the projection of the forcing into each 

634 eigenspace, gi(l) = ViH(l) for each i and I. If gi(l) = for all I, there is no 
ess fluctuation in the subspace containing eigenvector u,, that is, no fluctuation in 
636 Wi. In general, the more H is aligned with certain eigenvectors Uj, the more the 
es? fluctuations caused by the forcing signal will be concentrated in those subsys- 

638 terns. 

639 Similarly, a particular measurement generally observes some subsystems 

640 more than others. Assume we are observing the population via a scalar mca- 

641 surcmcnt, whether annual total population, recruitment or catch, represented 

642 as M(t) = Q(x(t)). Let us assume that Q(x) is differentiablc at the fixed point 

643 x. In the weak-noise limit this quantity also can be described by a linear ap- 

644 proximation, 

Q(x(tj) = + E f^) + °(M 2 ) 

1 Xl (A.6) 
= Q(x-) + J2<liyi(t) + 0(\y\ 2 ). 

i 

645 The measurement M(t) is approximated by the linear quantity Q{x)+J2i Qi Vi(t), 

646 which differs by only a constant from J2i Qi x i (*) i an d fluctuation in either quan- 

647 tity has the same characteristics as fluctuation in Y^iQiVi{t)- We can repre- 

648 sent the linearized measurement as a vector product qy(t), using a row vector 
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q = (qi, . . • , q n )- Changing to the natural coordinates of J, 

qy(t) = q^Wijt) Uj = ^qUiWi(t) = ^m t w t (t). (A.7) 

i i i 

Then rrii = qui expresses how strongly the motion in the subsystem containing 
Wi is reflected in the measurement. 

Frequency Analysis 

For the frequency analysis, take w(t) = and = for all t < 0. Then 



we may write the Z-transform (jElavdi 

oo 

w(z) = 53 w(k) z~ k 



1999) 



fe=0 



Atf(fc-i)+£(?(0£(*-o 

k=0 \ 1=0 / 

oo L oo 

= A ]T w(k 1) z- fe + £ G(0 £ £(* - I) ^- fc . 



(A- 



k=0 



1=0 k=0 



The bracketed expressions are themselves Z-transforms of shifted seque nces of 
w and £, so that by the shift property of the ^-transform (lElavdil . 11999). 

L 



){z) = A z- 1 d>(2) + 5] G(0 Z - 'f («) 



(A.9) 



657 or 



S ) = (i-ao _1 E^(o^» 



Z=0 

(1 — Az -1 ) -1 is a diagonal matrix with entries 
the ^-transformed deviation vector, 

L 

y(z) = Uii(z) = 53 



(A.10) 
r = Consequently 



z — A, 



is a linear sum of rational functions of z, with peaks tending to be near the 
eigenvalues of J. 
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662 Similarly, in the frequency analysis of a measurement M(t), 

oo 

M(z) =^M{k)z- k 

k=0 

K^{Q{^)+^m l w l {k))z- k (A.12) 

k i 

k i 

663 It is appropriate to discard the constant part of M(t) since we are concerned 

664 with year to year variation: 

i i 

= T M (z)i(z). 

665 The frequency response Tm{z) of the measurement M(t) is a weighted sum of 

666 the frequency responses for each subsystem, Ti(z) — jzr\., each weighted both 
ee? by the strength of environmental forcing in that subsystem at each lag, gi{l), 
668 and by the "emphasis" of the subsystem in the measurement, mj. 



L 
1=0 



z-Xi 



(A.13) 
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